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This talk is dedicated to the memory of Paco Yndurain, the original speaker in the conference. After a 
short account of his scientific career, we briefly review our ongoing collaboration to determine precisely the tttt 
scattering amplitude including the most recent data by means of Forward Dispersion Relations and Roy Equations. 
A remarkable improvement in precision over the intermediate energy region is obtained by using once-subtracted 
Roy Equations in addition to the standard twice-subtracted ones. 



1. Paco Yndurain, in memoriam 

F. J. Yndurain, "Paco" for his friends, was very 
dear for this QCD96-QCD08 Conference series, 
since, apart from an active participant, he was 
member of the Advisory Committee in all edi- 
tions. For the present QCD08, he was to review 
our work, but, very unfortunately, he passed away 
just before the beginning of the Conference. We 
thank prof. S. Narison for offering us the chance 
to give the same talk together with a brief account 
of Prof. Yndurain outstanding scientific career. 

In the University of Zaragoza, Spain, Paco ob- 
tained his Master degree in Mathematics in 1962 
and defended in 1964 his PhD thesis in Physics, 
"Definitions of Hamiltonians and Rcnormaliza- 
tion" , supervised by A. Galindo. In Zaragoza 
he also held an assistant professorship (1964/66), 
before moving to New York University as a Ful- 
bright Fellow (1966/67) and Associate Researcher 
(1967/68). He then became Research Fellow at 
CERN (1968/70), an institution to which he was 
specially attached, being senior scientific Asso- 
ciate in 1976 and 1985, and elected Member of the 
Scientific Policy Committee from 1988 to 1994. 

Finally, after a short period in the Complutense 
University of Madrid in 1970, he became full 
professor in the Autonoma Univesity of Madrid, 



where he had a decisive influence in shaping one 
of the best Theoretical Physics Departments and 
Particle Physics Groups in Spain. Actually, he 
was Director of the Department twice (1974/77- 
1981/84), dean of the Science Faculty in 1975 and 
research vice-rector (1978/81). 

Paco's intense activity was rewarded with 
many friends and collaborations all over the 
world, becoming visiting researcher at CERN, 
Michigan, Marseille, NIKHEF, Brookhaven, Or- 
say, Saclay, Viena, La Plata, Bogota, Caracas... 
He was also founding member of the European 
Physical Society, invited member of the Ameri- 
can Association for the Advancement of Science, 
the New York Academy of Science and both the 
Mathematical and Physics Spanish Royal Soci- 
eties. In addition, he became elected member 
of the European Board of High Energy Physics 
of the European Physical Society (1983/89), and 
elected president of the Physics and Chemistry 
section of the Spanish Royal Academy in 2002. 
But his influence extended beyond Physics as sci- 
entific advisor of IBM (1983/85), of the Institute 
for Scientific Research of Kuwait in (1980/82), 
and as elected member, by the Spanish Senate, of 
the Spanish University Council, as well as elected 
member of the World Innovation Foundation in 
London (2001) and the Academia Europaea. 
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Paco was awarded the gold medal of the Span- 
ish Royal Academy and was an honorific collab- 
orator of ICTP in Trieste, La Plata University, 
and Cavaliere Ufficiale nell'Ordine al Merito della 
Reppublica Italiana, a country very dear for him. 




1.1. Paco, Particle Physics and QCD 

Paco authored or coauthored more than 160 
publications, including articles and talks. Unfor- 
tunately in this space we can only highlight a few. 
In his New York period he was interested in an- 
alyticity methods for scattering theory, that he 
gathered in a Review of Modern Physics in 1972 
[T]. Around that time, QCD, the subject of this 
conference series, was established as the theory of 
strong interactions. This is where Paco obtained 
his most celebrated results, like the calculations 
of structure functions in the 70's and 80's j2], the 
hadronic contributions to the muon g-2 in the 80's 
and OO's |3j, quarkonium calculations in the 90's 
[4], as well as light quark and chiral symmetry 
QCD calculations in the early 80's [5]. 

The participants in this QCD08 conference 
might be more familiar with his QCD work, but 
he also made other well known contributions to 
the Standard Model phenomenology, like the cal- 
culation of radiative corrections to WW scatter- 
ing [6]. He also made some well known phe- 
nomenology beyond the Standard Model [7] but 
always limited within his very cautious and criti- 
cal attitude towards Standard Model extensions. 



In recent years, motivated by the very recent 
and precise experimental results, he was look- 
ing back to some of his initial interests, about 
analytic properties of pion-pion scattering ampli- 
tudes [819] ■ the subject of this talk. 

1.2. Books 

Paco dedicated an enormous effort and enthu- 
siasm to communicate Physics and Science. Of 
special relevance for the participants in this QCD 
Conference is his famous book - a classic - "QCD: 
An introduction to the Theory of Quarks and 
Gluons" (Springer- Verlag 1983). For many of 
us in the young - well, maybe not so young - 
generations of Particle Physicists, Paco's book 
has been THE BOOK where to learn QCD. He 
kept improving it throughout four different edi- 
tions (from the second it is entitled "The The- 
ory of Quark and Gluon Interactions"). More 
recently he also wrote a book on "Relativistic 
Quantum Mechanics: and introduction to Field 
Theory" (Springer- Verlag 1996), and another one 
on "Mecanica Cuantica" (In Spanish, 2nd Ed. 
Ariel S.A. 2003). 

Less known to the international community - 
they are written in Spanish- are his books on 
Popular Science and "Speculation" , as he called 
them, whose translated titles are: "Unified Theo- 
ries and the constituents of matter" , "Who walks 
out there" , "Electrons neutrinos and quarks" and 
"The Challenges of Science" . 

1.3. Why the name "Paco" and farewell 

Anyone called "Francisco" in Spain can also be 
called "Paco". The reason is that St. Francesco 
de Assisi - San Francisco in Spanish-, was the 
founder of Franciscan monastic community and 
so he was the "Father of the Community" or "Pa- 
ter Comunalis" in latin. Abbreviated: Pa. Co. 

But in the case of F.J. Yndurain, "Paco" was 
much more than just a nickname. In dire times 
for Spanish Scientific research, he was able to pro- 
duce first level results and was able to gather a re- 
search group on Particle and Theoretical Physics 
that later on became the Theoretical Physics De- 
partment of the Autonoma University in Madrid, 
extending his influence and scientific attitude 
much beyond. With his career and example he 
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really became one of the "PAter COmunalis" of 
the Spanish Particle Physics Community. 

Paco, we will miss you as a friend, but also 
your enthusiastic and contagious vitality not only 
about Physics but life as a whole. Sit tibi terra 
levis — may the earth rest lightly on you. 

Of course, knowing Paco, I guess that by now 
he would be highly embarrassed and yelling from 
the back... Enough is enough! Let's have fun with 
Physics!!... - Let's not make him wait any longer. 

2. Dispersive approach to tttt scattering 

A precise determination of the tttt scattering 
amplitude at low energies is relevant for the study 
of Chiral Perturbation Theory (ChPT), quark 
masses, the chiral condensate [10_ and, at inter- 
mediate energies, for the properties of the contro- 
versial sigma meson. However, the existing exper- 
imental information from tttt scattering has many 
conflicting data sets at intermediate energies and 
for many years very little data in the low energy 
region close to threshold. The interest in the topic 
has been renewed with recent 11^ and precise ex- 
periments on kaon decays, relatively easy to relate 
to TTTT scattering. 

From the theoretical side, this process is very 
special due to the strong constraints from isospin, 
crossing and chiral symmetries, but mostly from 
analyticity. The latter allows for a very rigorous 
dispersive integral formalism that relates the tttt 
amplitude at any energy with an integral over the 
whole energy range, increasing the precision and 
providing information on the amplitude even at 
energies where data are poor, or in the complex 
plane. Remarkably, it is model independent since 
it makes the data parametrization irrelevant once 
it is included in the integral. The dispersive ap- 
proach is thus well suited to study the threshold 
region or poles in the complex plane associated 
to resonances. Our recent works [8 9 make use 
of two complementary dispersive approaches, in 
brief: 

• Forward Dispersion Relations (FDRs): Cal- 
culated at t = so that the amplitude unknown 
large-t behavior is not needed. There are two 
symmetric and one asymmetric isospin combina- 
tions, to cover the isospin basis. The symmetric 



ones, 7r'^7r+ and tt'^tt'^, have two subtractions 
YleF{s,0)-F{AMl,Q)= (1) 

s(s-4Af^) p p roo (2s'-4M^)Im_F(s',0)d^' 

where F stands for the Fo+(s,t) or Foo(s,t) am- 
plitudes. They are very precise since all the in- 
tegrand contributions are positive. The antisym- 
metric isospin combination /t = 1 reads 

i^(^*=i)(s,0) = (2) 

2s-4M|_pp r-x, , I ImF<^* = ^(a',0) 
TT J4M2 *^ (s'-s)(s'-|-s-4Af^) ■ 

AU FDRs are calculated up to - 1420 MeV 

• Twice subtracted Roy Equations (RE): they 
are an infinite set of coupled equations |13| . 
equivalent to nonforward dispersion relations plus 
some t — s crossing symmetry. They are well 
suited to study poles of resonances since they are 
written in terms of partial waves f^^^ of definite 
isospin I and angular momentum I. The compli- 
cated left cut contribution is rewritten as a series 
of integrals over the physical region: 

Refr\s) = c[%^^^+C[^%^^^ (3) 

+ Vp.R / ds%,z-a,7'(s',s)Im///'(s'), 

J, J AMI 

where the c'l^\ C';*'^'' constants and Ki^i-jji ker- 
nels are known. The calculation is truncated 
at / < 2 and at some cutoff energy sq. The 
/ > 2 waves and the high energy are input. Al- 
though RE are vahd for < 8M^ ~ 1120 MeV, 
for now we have implemented them only up to 
2Mk. 

The aim of our group when using RE com- 
bined with FDR has been to improve the pre- 
cision of scattering data analysis, to test ChPT, 
and, still in progress, to obtain precise determi- 
nations of the fj resonance. The Bern group [14j 
has also carried out a series of RE analysis, with 
and without ChPT constraints, using as input 
phenomenological parametrizations for the / > 2 
waves and above 800 MeV for the other waves, 
as well as some Regge input. When using ChPT 
constraints, they find, 4°^ = 0.220 ± 0.005 M^^ 
and = -0.0444 ± 0.0010 A/^\ an extremely 
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precise claim, together with predictions for other 
scattering lengths and the S and P wave phase 
shifts up to 800 MeV. Although of the original 
input, particularly the Regge theory and the D 
waves, was questionable [15j, it certainly seems 
to have a very small influence in the threshold 
region of the scalar waves [16j . In addition, the 
Krakow-Paris [17] and Paris [18] groups have per- 
formed other RE analysis. The former resolved a 
long-standing ambiguity, discarding the so-called 
"up" ' solution, including in their analysis a study 
using polarized target data. The latter checked 
the calculation in and claimed an small dis- 
crepancy in the Olsson sum rule. 

3. Sketch of our analysis 

The approach we have followed over a series of 
works [S] can be sketched as follows: 

(1) We obtain simple "Unconstrained Fits to 
Data" (UFD) to each tttt scattering wave so that 
it can be improved independently if needed. High 
energy fits use Regge theory. As our precision 
improves, we refine our fits with more flexible 
parametrizations or with new and more precise 
data, as, for example, in [9] where we included 
the newest and very precise Ki^ data [TT] . 

(2) We then check how well these UFD sat- 
isfy FDRs and RE and some sum rules. To 
our surprise some of the most widely used data 
parametrizations fail to satisfy FDRs or the sum 
rules. Thus we choose the data parametrizations 
in better agreement with FDRs. 

(3) Our best results are obtained by imposing 
dispersion relations (both FDRs and RE) and 
some crossing sum rules (SR) on the data fits, 
now called "'Constrained Fits to Data"' (CFD), 
which are thus consistent with analyticity, unitar- 
ity, crossing, etc... 

A dispersion relation i is well satisfied at a 
point Sn if the difference between the left and 
right sides in Eqs.([T]),([2]) and ([3]) is small relative 
to its uncertainty SAi. Thus, when the average 
discrepancy (similar to an averaged / {d.o.f)) 



* number of points \(5Ai(s„) J 

it implies that the corresponding dispersion rela- 
tion is satisfied within uncertainties. In practice, 



the values of s„ are taken at intervals of 25 MeV. 
This measure is also used [5] to obtain the Con- 
strained Fit to Data, by minimizing: 

/ cxp \ 2 

where p^^^ are all the parameters of the different 
UFD parametrization for each wave or Regge tra- 
jectory, thus ensuring the data description, and d/ 
and dj are the discrepancies for a couple of cross- 
ing sum rules. The weight W = ^ was estimated 
from the typical number of degrees of freedom 
needed to describe the shape of the amplitude. 

The result of this program is, on the one hand, 
a set of precise Constrained Fits to Data, that 
satisfy very well all dispersion relations within un- 
certainties. Remarkably, all the waves are given 
in terms of very simple parametrizations that can 
be found in [8] which are very easy to use for phe- 
nomcnological purposes. On the other hand, we 
have the outcome of Roy Eqs. themselves that al- 
lows us to extend the calculation to the complex 
plane and look for poles associated to resonances 
and study their parameters. 

In particular the best determination of thresh- 
old parameters is obtained by using the CFD set 
directly or inside appropriate sum rules 8]. For 
the SO and S2 waves, we find: a\^^ = 0.223 ± 
0.009 and = -0.0444 ± 0.0045 M"!, 

in remarkable agreement with the predictions in 
[T^ . that extends also to the P wave scattering 
length. There are however, some disagreements 
of 2 to 3 standard deviations in the P-wave slope, 
and also in some D wave parameters. In general, 
for the controversial SO wave, the agreement is 
fairly good only up to roughly 450 MeV, but from 
that energy up to 800 MeV those predictions de- 
viate slightly from our data analysis. This means 
that we find from our Roy Eqs. and FDR data 
analysis a somewhat different pole from that of 
the Bern group (see R. Kaminski talk in this con- 
ference). Let us emphasize that we are talking 
about a deviation of a few degrees that only af- 
fects the sigma mass and width determination by 
ten or twenty MeV at most, which is a remarkable 
improvement compared with the situation just a 
few years ago and the huge and extremely conser- 
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Figure 1. Above 500 MeV, the uncertainty in stan- 
dard SO Roy Eqs. (Top) is much larger than with 
once subtracted Roy Eqs. (Bottom). 

vative uncertainties of hundreds of MeV quoted 
in the PDG at present pj] . 

3.1. Latest developments 

In this conference we report prehminary re- 
sults on two issues: First, we have recently im- 
plemented once-subtracted Roy Eqs, (denoted 
GKPY for brevity) [20]. This is motivated by the 
interest in the 450MeV < y/s < 2Mk region 
of the SO wave, dominated by the cr-resonance. 
The large low-energy experimental uncertainties 
in the S2-wave translate into large uncertainties 
for the dispersive output for the SO wave in that 
region when using the standard twice subtracted 
Roy Eqs. (shaded area in top panel of Fig 1), but 
not for the once subtracted ones, which have a 
much smaller uncertainty (Fig 1, Bottom). 

Second, we have improved the matching be- 
tween the SO wave parametrizations at low and 
intermediate energies, which occurs at 932 MeV, 
imposing also continuity in the first derivative. 




s"' (MeV) 




400 600 800 1000 

s""(MeV) 



Figure 2. Top: The improved matching only changes 
slightly the SO wave above 932 MeV. Bottom: The 
agreement GKPY Eqs. is slightly improved above 
850 MeV (compare with Fig.l Bottom). 

This follows a suggestion ^21j to explain the 
roughly 2a level discrepancies in the SO wave be- 
tween our analysis and that of the Bern group in 
the 400 — 900MeV region. In Fig. 2 we compare 
the CFD in [8] (dashed line) and the new one with 
the improved matching (continuous line), which 
only differ a little above 932 MeV. Note that, for 
clarity, we do not provide data points, which are 
nevertheless reasonably well described by both 
parametrizations when taking into account exper- 
imental errors. Hence, the disagreement is not 
caused by a poor matching, but we are studying 
further suggestions in [51] . Nevertheless this bet- 
ter matching improves the Roy and GKPY Eqs. 
fulfillment above 850 MeV, as seen by comparing 
the bottom panels in Figs.l and 2. 

Finally, in order to show how well our prelim- 
inary CFD set satisfies the nine dispersion rela- 
tions, we list in Table 1 the average discrepancies, 
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FDRs 


< 932 MeV 


s^^^ < 1420 MeV 




0.13 


0.31 




0.83 


0.85 


It=i 


0.13 


0.70 




< 992 MeV 


Roy Eqs. 


2 subs. 


1 sub. (GKPY) 


SO 


0.06 


0.49 


S2 


0.13 


0.11 


P 


0.11 


0.23 



Table 1. Average discrepancies of the Constrained 
Fit to Data for each forward dispersion relation. 



defined in Eq.Q, for the three FDR, up to two 
different energies, and for the three Roy Eqs. up 
to ^ 2Mk, either with the standard two subtrac- 
tions or only with one. From that Table it is 
clear that this CFD set satisfies remarkably well 
all dispersion relations within uncertainties 

In summary the CFD set provides a model in- 
dependent and very precise description of the tttt 
scattering data consistent with analyticity and 
crossing that can be easily used for phenomeno- 
logical purposes. Work is in progress to obtain an 
even more accurate description for threshold pa- 
rameters than the one we had, and (see R. Kamin- 
ski talk) by using it inside once-subtracted Roy 
Eqs., to obtain a precise determination from data 
of the cr pole parameters. 
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